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Spin squeezing of a nonlinear interaction model with Josephson-like coupling is studied to obtain 
time scale of maximal squeezing. Based upon two exactly solvable cases for two and three particles, 
we find that the maximal-squeezing time depends on the level spacing between the ground state 
and its next neighbor eigenstate. 
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I. INTRODUCTION 



The phenomenon of spin squeezing in collective spin 
system have attracted much attention for decades not 
only because of fundamental physical interests 0, 0, H, 
EL B IE S IE @) but also for its possible application in 
atomic clocks for reducing quantum noise [E IE 0, H| and 
quantum information [IE [Tl|, [IE EE Hi- The occur- 
rence of spin squeezing is due to quantum correlations 
among individual spins, which requires at least two spins 
and nonlinear interaction between them. Kitagawa and 
Uea have studied the spin squeezing generated by the 
so-called one-axis twisting (OAT) model with Hamilto- 
nian: Hoat = 2k III- Possible realization of the OAT- 
type squeezing in a two-component Bosc-Einstein Con- 
densate (TBEC) (lOL |l 5l| , and atomic ensemble system in 
a dispersive regime [16| have been investigated recently. 
S0rensen et al. also proposed that the spin squeezing can 
be used as a measure of many-particle quantum entan- 
glement [Tol | . 

So far the OAT-type spin squeezing was mainly studied 
in Heisenberg picture. As a result, the explicit expres- 
sion of the spin squeezed state is unknown. Moreover, 
the direction that spin squeezing is observed varies with 
time [ll|. Jaksch et al. have shown that the OAT-type 
SSS can be stored for arbitrarily long time by removing 
the self- interaction [l7|. However, it might not be easy 
to handle in experiment since the precisely desi gned ad- 
ditional pulses are crucially required. In Refs. fill fTsj] . 
the authors proposed the constant-coupling scheme by in- 
troducing additional Josephson-like coupling fl J x to the 
OAT model. It was shown that the Josephson interaction 
results in an enhancement of spin squeezing compared 
with that of the OAT. Moreover, the strongest squeezing 
appears in the z direction [l8j |. which, however, is valid 
only at the maximal-squeezing time (MST). Although 
some formulas of the MST for extremely small [I| or large 
coupling [IE [H, HE [U, HE] have already been known, it 
is challenging to determine the MST within an interme- 
diate coupling I < 0/ k << N [where N is total particle 
number] . 

In this paper, we reconsider the constant-coupling 
scheme [IE HE] with the purpose to determine the MST. 
We find all the analytic solutions for two- and three- 



particle cases. Motivated by the exactly solvable cases, 
we show that the MST depends on the level spacing be- 
tween the ground state and its next neighbor eigenstate. 
We explain it by investigating the spectral distribution of 
the spin state, and find only the two lowest available lev- 
els are predominantly occupied. Our paper is organized 
as follows. In Sec. II, we introduce theoretical model and 
derive some basic formulas. To proceed, in Sec. Ill, we 
gives some analytic expressions for the cases of N = 2 
and N = 3. In Sec. IV, we study the spin squeezing 
for many-particle cases, and present exact diagonaliza- 
tion method to obtain the MST. Moreover, we compare 
our result with its analytic solution. Finally, a summary 
of our paper is presented. 



II. THEORETICAL MODEL 



Formally, a two-level atom can be regarded as a 
fictitious spin- 1/2 particle with spin operators — 

(\b)u(b\ - \a) u (a\)/2, and sf = (a<?)t = |6>«<o|, where 
\a)i and \b)i are the internal states of the ith atom. We 
consider an ensemble of N atoms with its dynamics de- 
scribed by collective spin operator: J = J2t-i s ^ ■ T nc 
spin squeezing is quantified by a parameter [if: 



V2(AJ n ) 



J 



1/2 



(1) 



where j = N/2, and (A j n ) m in represents the smallest 
variance of a spin component J n = J • n normal to the 
mean spin (J). For a coherent spin state (CSS), the 
variance (AJ n ) m j n = \J jjl and £ = 1. In general, a spin 
state is called spin squeezed state (SSS) if the variance of 
the spin component J n is smaller than that of the CSS, 
i.e. i < 1. 

Follow Refs. (H Hi HI SI, M, El H [3(3 , we consider 
a nonlinear spin system governed by 



h = nj x + 2kJ; 



(2) 



which can be realized in the TBEC [3J, |3J. The first 
term is Josephson-like coupling induced by a microwave 
(radio frequency) field. The Rabi frequency f2 can be 



2 



controlled by the strength of the external field. The sec- 
ond term is the self-interaction aroused from nonlinear 
collision between atoms. An initial coherent spin state 
— e~ l7rJy / 2 \ j, — j) will be considered in this pa- 
per. Physically, the Dicke state \ j, —j) represents all the 
atoms occupying in the internal ground state |a). By 
applying a short ir/2 pulse to the Dicke state, one can 
obtain the CSS with each spin to be aligned along the 
negative x direction 10]. After that, one switches on the 
Josephson-like immediately, then dynamics of the spin 
system is governed by the Hamiltonian Note that, we 
will consider only positive k case. However, our results 
keep valid in the opposite case by using initial maximum 
weight state of J x , i.e., \j,j) x - 

The state vector at any time t can be expanded in 
terms of eigenstates of J z : \ip{t)) — ^ m c m (t) \j,m) , 
where — j < m < j. The probability amplitudes c m (t) 
can be solved by time-dependent Schrodinger equation, 
obeying 

— E 7n C 7n -\- X m C rn — i -\- A_ m C m _)_i, (3) 

2«;m 2 , and X m 

A 3 



where E. 
with A_ 



— 



+ m)(j -m + 1) 
j — and X± m = X^ m+ \. The probability 
amplitudes of the initial CSS 



,(0) 



(-1) 



j+m 



(2j) 



2J 



(j - m)\(j + m)!' 



(4) 



satisfy c_ m (0) = c m (0) for even N, and c_ m (0) = 
— c m (0) for odd N. Due to the symmetry properties of 
the elements X± m and the initial amplitudes c m (0), we 
obtain simple expressions: C- m (t) = ±c m (t), which in 
turn result in (J y ) = (J z ) = 0, and (J x ) ^ 0, i.e., the 
mean spin ( J) is always along the x axis. The spin com- 
ponent normal to the mean spin is J n = J y sin 9 + J z cos 9 
and its variance is (AJ n ) 2 = (J 2 ) — (J n ) 2 = \C + 
co|2e A+ sir^2£ Bj where A = (J 2 — J 2 ), B = (J z J y +J y J z ), 

and C — ( J 2 + Jy). By minimizing the variance (A J n ) 2 
with respect to 9, we get the squeezing angle: 

1 



•tan _1 (B/A), 



(5) 



and the smallest variance 

(AJ^^C-iv^T^, (6) 

from which one also obtain the squeezing parameter Eq. 
(frj. We consider the spin squeezing in the intermediate 
coupling regime, namely 1 < Q/k << N, where no an- 
alytic solutions are available for the nonlinear spin sys- 
tem [Til f33j | . However, we can exactly solve two- and 
three-particle cases. Some of important physics can be 
extended to many-particle cases. 



III. EXACT SOLVABLE CASES 

In this section, we study the spin squeezing based on 
two exact solvable cases with N = 2 and N = 3. Though 



simple, it is of general interest to investigate the relation- 
ship betwee n sp in squeezin g an d quantum entanglement 
03 El El, El H SE [M ll- Such a relationship for 
two-particle (^two-qubit) |3l. FIT 35L l36j and three-particle 
[U, [33| have been studied recently. Here, we focus on 
dynamical behavior of the spin system to show the con- 
ditions of the optimal squeezing and its time scale. 



A. Two-particle case 

For the simplest case N = 2 (j = 1), only three spin 
projections (m = —1, 0, +1) are involved. From Eq. ([3]), 
we obtain 




)■ 
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2X t \ 
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(7) 



where E m and X m are defined in Eq. ([3]) , and we have in- 
troduced the linear combinations of the probability am- 
plitudes p[ + \t) = d{t) + c_i(t) and p^\t) = 2c (t), 
with the initial conditions p^ + ^(0) = 1 and Pq^\o) = 
—y/2. Similarly, we also introduce Pi(t) — ci(t) 
i(t). However, its solution p^ - ' (t) 



c-i 



-i2«t„(— ) 1 



due to p[ \0) = 0. Therefore, we obtain c\{t) = 
c_i(t) = p[ +) /2, which gives (J z ) = |c +1 | 2 - |c_i| 2 = 
and (J + ) = V2(c_ic5 + c c|) = 2\/2Re(coc 1 ). Since (J + ) 
is a real function, ( J y ) = and (J x ) ^ 0, which show 
that the mean spin is always along the x direction. Such 
a result is valid for arbitrary even N . Eq. jJJ can be 
solved exactly, then one obtain immediately the reduced 
variance 



e-"~™ -p\ '(0) = 



(A J n )Ln = \ ~ I Ism St\ ^1 - ^ sin 2 St, (8) 



and the squeezing angle 



'min — 2 ^ an 



Scos(St) 
tt sm(St) 



(9) 



where 25 = £ 3 — £ 1 = 2\/^ 2 + « 2 is the level spacing 
between the second excited state £3 and the ground state 
£\ , obtained by solving the eigenvalues of the coefficient 
matrix of Eq. ((7]). As shown in Fig. [T|a), we find that at 
the times t\ = kn/S, £ revives periodically to its initial 
value 1. In fact, apart from a globe phase, the states 
at i*, \ip(t* k )) = {-l) k e - iKt *x 11,-1)3., are just the initial 
CSS. 

From Eq. ©, we find that the vanishing 9 m [ n occurs 
at tk = (k+ l/2)7r/5, and the state vector at tk reads 

|V(t*)) = {-Ifie' 1 ^ 

xjsin^lE-^^-cos^lE+^J^lO) 

which correspond to a superposition of two coherent spin 
states, (1,-1)2, and |1, +l)a; with the mixing angle 77 = 



3 



tan (Q/k). Obviously, if the coupling is very strong 
(O ^> k), sin(^) = n/S — ► 1 and cos(i]) = k/S — > 0, so 
IV'(^fc)) ~~ * |1, — which in turn leads to a very weak 
squeezing at tk- On the other hand, if the coupling is 
very weak (fi <C n), \ip(tk)) — > |1, l)a:, which also results 
in a weak squeezing at tk- Therefore, we will study the 
spin squeezing within the intermediate coupling regime. 



1.0 



5 0.5 



0.0 

it/4 



min 0.0 



-Jt/4 L 

0.0 





„ " "* N 

Q=5k 




< a k - 

o_=2k 

fi=K 


\ 

>. 












... " 

\ 1 " - . 
\ 1 




\ 

s 

\ 






S 




\ | 

1 



0.5 



1.0 

Kt 



1.5 



2.0 



FIG. 1: (color online) Time evolution of (a) the squeezing 
parameter and (b) the squeezing angle for N = 2 and various 
Rabi frequencies: fl = 5n (Dashed red lines), fl = 2k ( Dotted 
blue lines), and Q = k (Solid black lines). The maximal- 
squeezing times to for different Q are indicated by the vertical 
lines. 



B. Three-particle case 



For N = 3 ( j = 3/2) case, we introduce the linear 
combinations of the amplitudes pffl (i) = c m (t) + c_ m (t) 



with m = 1/2,3/2. Since p[V 2 (0) = Pul(°) = °> we S et 



,(+)< 



3/2 » 



'1/2 * 



Pm" (t) = 0. Therefore, the amplitudes obey c m (t) — 
— c- m (t), from which we can prove that the mean spin 
is always along the x direction. Such result keeps valid 
for any odd N case. From Eq. ([3]), we obtain a coupled 

equations for the linear combinations pL \t) — c m (f) — 
c- m (t): 



(12) 




where E[, 2 = E 1/2 - X 1/2 - 
,(-)/ 



The initial conditions are 



Pyl(6) = — 1/\/2 and py^O) = Dynamical evo- 

lution of the three-spin system is determined solely by 
Eq. (fT2| . The analytic expression of the variance is 



(AJ n ) 



2 

min 



3k 



3 

4 + 52 
3 k 



sin 2 St 
sin 5i| 



1 



3k 2 
S 2 " 



sin 2 St, (13) 



and the squeezing angle is 
1 



'min — 2 ^ n 



ScoajSt) 
(n + k) sm(St) 



(14) 



In Fig. [TJ time evolutions of £ and min are investigated 
for the coupling f2 > k. We observe that local minima of 
£ together with 6* m i n = also occurs periodically at the 
times tk- Moreover, with the decrease of f2, the squeezing 
parameter at tk becomes small, i.e., more squeezed. For 
the coupling f2 = n, the spin system is optimally squeezed 
at tk , as shown by the solid black lines of Fig. [TJ In this 
case sin(7y) = cos(?7) = l/v2, and the spin states at tk 
are 



m k )) = (-1) 



. ie 



{11,-1) 



V2 

i(-l) k+1 e- iKtk \j = l,m=0). (11) 



Here, the state (|1, — l) x — \l,+l) x )/y2 is maximally 
entangled (or Bell) state, while the Dicke state \j = 
1 , m = 0) is maximally squeezed state [2j . For this 
state, both the mean spin (J x ) and the variance (AJ n ) m j n 
are equal to zero, which makes it hard to define £ as 
Eq. UJ). To avoid this problem, Wineland et al. proposed 
another definition of the squeezing parameter, namely 
£ — > C// 1 (</)!)£, which gives the smallest squeezing l/v2 
for TV = 2 case 0. 



where 2S = £3 - £1 = + 2kH + 4k 2 is the level 

spacing for N — 3 case, and £3 and £\ correspond to two 
eigenvalues of the coefficient matrix of Eq. ([12")) . 
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FIG. 2: (Color online) Time evolution of (a) the squeezing 
parameter, and (b) the squeezing angle for N = 3 case with 
various Rabi frequencies: Q = 5k (Dashed red lines), Q — 2k 
(Solid black lines), and fl = k (Dotted blue lines). 

In Fig. [51 we investigate time evolution of £ and # m i n 
for N — 3 case. Similar with previous N = 2 case, our 
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results show that for 51 > 2k, local minima of £ together 
with 6* m i n = also occur at the times S% = (k + l/2)n. 
As shown by the solid black line of Fig. [2j we find that the 
optimal squeezing can be obtained at tk for the coupling 
51 = 2k. The maximally squeezed state at tk reads 



m k )) = 



»(-!)* 
V2 



3i«t fc /2 ^ 
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0- 
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1 




2' 




2' 


2 



(15) 



Such a state gives the smallest squeezing parameter 
£,{tk) — 1/V3 that the three-particle system can reach. It 
is worth mentioning that for 51 < 2k, the vanishing m i n 
appearing at tk no longer corresponds to local minima 
of £, as shown by the dotted blue lines of Fig. O The 
time scale tk is relevant to determine the MST only for 
51 equal or larger than the optimal coupling. 

In short, we find some basic features for two exactly 
solvable cases. Local minima of £ with # m i n = occur 
at the MST tk- This is no longer true if 51 smaller than 
the optimal coupling. The time scale tk depends on the 
level spacing 2S between the ground state and the second 
excited state. Due to the symmetric properties of the spin 
system, the first excited eigenstate is an idle level (see 
below) . This is also the reason why we can introduce the 

linear combinations of the amplitudes Pm , with to — 
0, 1 for N = 2, and m = 1/2,3/2 for N = 3. For the 
optimal coupling, the spin system will be evolved into 
the maximally squeezed state at tk ■ 1 1 , 0) (for N = 2) 
or (|3/2, 1/2) - |3/2, — l/2>)/v / 2 (for N = 3), which is 
just the ground state of 2k J%. We will extend the above 
results to many-particle cases. 



IV. MANY-PARTICLE CASES: THE 
MAXIMAL-SQUEEZING TIME 

In this section, we study the spin squeezing for many- 
particle cases focusing on the time scale of the maximal 
squeezing. For instance, we consider the spin system with 
particle number TV = 40 [j], [l6| ■ The numerical results 
are shown in Fig. [3j We find that with the increase of 51, 
the squeezing £ and the mean spin ( J x ) show collapsed 
oscillations [H, [38|. Local maxima of the mean spin (J x ) 
always appear together with the vanishing # m i n . We can 
prove this from Heisenberg equation of J x and Eq. ([5]): 



d(J x )/dt 



(JzJii 



JyJz 



A tan(20 m i„). If the mean 



spin reaches its local maximum at a certain time to, then 



d(J x )/dt 



0, which leads to 9 n 



at to provided 



that A f 0. 

As shown in Fig.[3jb), for a small coupling with 51 = k, 
there are two time scales: to for the vanishing # m i n , and 
To for the maximal squeezing. Note that the latter time 
scale r closes to that of the OAT result (51 = case), 
i.e. ktq ~ 0.04986. With the increase of Q, these two 
time scales become coincident, as shown in Fig. ^c) and 
(d). Unlike to the exactly solvable cases, we find that 
the optimal coupling for TV = 40 is not a fixed value but 
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FIG. 3: (color online) Time evolution of £ (thick solid lines), 
#min (dashed blue lines) and (Jx) / j (thin red lines) for N = 40 
(j — 20) and various Rabi frequencies: (a) 51 = 0, (b) 51 = k, 
(c) 51 = 4.24k (the optimal coupling), and (d) 51 = 20k. The 
time scale to for different 51 are indicated by the vertical lines. 



can be arbitrary 51 in a region 4.239 < Qr/k < 4.242. 
Fig. Etc) represents the optimal squeezing case with the 
coupling 51 = 4.24k. Starting from the initial CSS, the 
spin system evolves into the maximally squeezed state at 
Kt = 0.09. 

To investigate the maximally SSS at to, we calculate 
the quasiprobability distribution (QPD, or the Husimi 
function) Q(6,(f) on the Bloch sphere [l| 



Q(M)H(0,</#(i)>| 2 , 



(16) 



where \0,<j>) = exp{— i0( J x sm<j> — J,, cos <j>)}\i,—j) is the 
generalized coherent spin state [3^, |4(| • The initial state 
is a particular case of the CSS, namely \j, —j) x = \9 = 
ir/2,<fi = 7r). The QPD can be used to simulate the 
variation of spin uncertainties. The circle in Fig. H^a) 
represents an isotropic spin variance for the initial CSS, 
while the shaded ellipse parts in Fig. 0Jb) an d (c) are 
that of the SSS at times about to/2 and to, respectively. 
Unlike to the OAT result the maximal variance 

reduction appears along the z axis with 6* m ; n = [l8l ]. 
In Fig. [H we also calculate the probability distribution 
\c m \ 2 — m\ip(t))\ 2 of the spin state for N — 40 and 
the optimal coupling 51 = 4.24k. Compared with the ini- 
tial CSS, we find that the maximally SSS at to has a very 
sharp probability distribution with a large amplitude of 
the lowest spin projection, i.e., m = (for even N) or 
to = ±1/2 (for odd TV) [4l|. Such a sharp probability 
distribution of the SSS can be explained qualitatively by 
considering the familiar phase model [42j (see also refer- 
ences therein). 

In order to determine the time scale to, we employ nu- 
merical diagonalization of the Hamiltonian @ to obtain 
a set of eigenenergies {S n ;n = 1,2,3,...}, where n = 1 
denotes the ground state, n = 2 the first excited state, 
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FIG. 4: Time evolution of the quasiprobability distribution 
Q(8,(f)) (top) and the probability distribution jc m | 2 (bottom) 
at the times (a) nt = 0, (b) nt = 0.04, and (c) nt = 0.09 (the 
MST). The QPD is normalized such that Q(tt/2, tt) = 1. In 
(b), the spin component normal to the x axis is defined as 
J n = J • n with n = (0, sin 0, cos 8). Other parameters are 
taken as those of Fig. [3jc) . 

and n = 3 the second excited state, etc., thus 

H\ct> n }=£ n \4> n ), (17) 

where £ n and \(f> n ) with n — 1, 2, 3, • • • , (2j + 1) depend 
on the parameters N and Q [43j |. As shown in Fig. [5fa), 
we plot parts of £ n for j = 20 (N = 40) as a function 
of the coupling f2. Similar with previous two- and three- 
particle cases, we suppose that the MST t depends on 
the level spacing between n = 1 and n — 3, namely 
t = 71-/(25) with 2S = £ 3 - £ x . To check it, in Tablell 
we compare exactly numerical results of the time to with 
7i"/ (25) for various N and fi. Our diagonalization method 
gives accurate prediction of the time to- We remark that 
for N — 2 and N — 3 cases, both two results are exactly 
the same. 

To explain the above agreements, we calculate the 
spectral distribution of the spin state, i.e., \((f> n \ip(t))\ 2 in 
Fig. OHb)-(d) for N — 40 and various f2. Physically, the 
spectral distribution measures the population distribu- 
tion of the state vector \ip(t)} on the eigenstates \4> n ) 0- 
For fixed parameters N and fl, the spectral distribution 
| (4> n \ip(t)} | 2 is time-independent. In fact, one can expand 
the spin state in terms of {|0 ra )}: \?P{t)) = J2 n d n {t)\4>n) 
with the amplitudes d n (t) = exp[— i£ n t]d n (0). Here 
the initial amplitudes d n (0) depend only on the ini- 
tial condition Eq. (|4|), therefore the spectral distribu- 
tion |(0„|VW)| 2 = K(t)\ 2 = K(0)| 2 and is time- 
independent for fixed N and Q. From our numerical 
calculations, Fig. E]Jb-d) , we find that total occupation of 
the spin state \ip(t)) on the eigenstates n — 1 and n = 3 is 
over 80 percent. This is the reason why the MST depends 
on the level spacing between these two levels. Moreover, 
we find the even n eigenstates are in fact the idle levels, 



just as previous N = 2 and N — 3 cases. 

Except for N = 2 and N = 3, exact solutions for the 
nonlinear spin system within the small- coupling regime 
(1 < tt/n « N) do not exist [H, HI]. In our previ- 
ous work |4lj |. however, we have obtained the analytic 
expression of the MST based upon the phase model: 

7T / K 

^Waw' (18) 

which is valid for large N (> 10 3 ). Our analytic solu- 
tion of the MST is derived by the prediction to ~ T/4, 
where T = 2tt /ui e s is the period of the pendulum near the 
bottom of a periodic potential [4l| . In fact, for large N 
the spin system behaviors as a pendulum rotating with 
the oscillating frequency cj c ff = ^/2kQ,rN . As shown in 
Table HI we compare ir/(2S), the analytic solutions of 
Eq. (US]), and the exact numerical results of the MST 
for various parameters f2 and N. It is shown that our 
analytical expression of Eq. (|18p works very well for the 
large TV (~ 10 3 ), which implies that the oscillating fre- 
quency w c ff has its physical meaning to be half of the level 
spacing S — (£3 — £i)/{2K). Note that the phase model 
or Eq. (fT5|) is valid for the large N, while w/(2S) is no 
limited by this. From this sense, we believe that the diag- 
onalization method presented here provides much com- 
prehensive way to measure the maximal-squeezing time. 




FIG. 5: (color online) (a) Part of the Eigenenergies £ n as 
a function of Q for N = 40. The spectral distribution 
\((f>n\ip(t))\ 2 for (b) fl = k, (c) fl = 4.24k, and (d) SI = 20k. 



V. CONCLUSIONS 

In summary, we have studied the maximal-squeezing 
time of a nonlinear spin system, which can be realized 
in the two-component BEC, or other spin system similar 
with Takeuchi et al.,16]. Motivated by two exactly solv- 
able cases for N = 2 and N — 3, we show that time scale 



6 



TABLE I: Comparison of exactly numerical to, tv/(2S), and analytic results of Eq. (I18|l for different N and SI. The times are 
in units of (100k) _1 . 





N = 40 


N = 200 


N = 1000 


Q/k: 


1 4.24 20 


1 6.7 25 


1 10.8 50 


Exact num. 
n/[2S) 
Eq. (18) 


18.28 9.065 3.604 
19.02 8.615 3.573 
17.56 8.529 3.927 


8.192 3.184 1.549 
8.143 3.048 1.533 
7.854 3.034 1.571 


3.665 1.104 0.4945 
3.571 1.071 0.4916 
3.512 1.069 0.4967 



of the maximal squeezing depends on the level spacing 
between n = 1 and n = 3 eigenstates. We explain it by 
calculating the probability distribution of the spin state 
on the eigenstates of the Hamiltonian, and find that the 
above two states are occupied predominantly. Such re- 
sults keep valid for arbitrary N and a wide rage of the 
coupling strength. 
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